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Abstract 

Two series of uniaxial tensile tests are performed on isotactic polypropylene with 
the strain rates ranging from 5 to 200 mm/min. In the first series, injection- molded 
specimens are used without thermal pre-treatment. In the other series of experiments, 
the samples are annealed for 51 hour at 160 °C prior to testing. 

A constitutive model is developed for the viscoplastic behavior of isotactic polypro- 
pylene at finite strains. A semicrystalline polymer is treated as an equivalent hetero- 
geneous network of chains bridged by permanent junctions (physical cross-links and 
entanglements). The network is thought of as an ensemble of meso-regions connected 
with each other by links (lamellar blocks). In the sub-yield region of deformations, 
junctions between chains in meso-domains slide with respect to their reference posi- 
tions (which reflects sliding of nodes in the amorphous phase and fine slip of lamellar 
blocks). Above the yield point, the sliding process is accompanied by displacements of 
meso-domains in the ensemble with respect to each other (which reflects coarse slip and 
fragmentation of lamellar blocks). To account for alignment of disintegrated lamellar 
blocks along the direction of maximal stresses (which is observed as strain-hardening of 
specimens in the post-yield region of deformations) elastic moduli are assumed to de- 
pend on the principal invariants of the right Cauchy-Green tensor for the viscoplastic 
flow. 

Stress-strain relations for a semicrystalline polymer are derived by using the laws 
of thermodynamics. The constitutive equations are determined by 5 adjustable pa- 
rameters that are found by matching observations. Fair agreement is demonstrated 
between the experimental data and the results of numerical simulation. A noticeable 
difference is revealed between the mechanical responses of non-annealed and annealed 
specimens: (i) necking of samples not subjected to thermal treatment precedes coarse 
slip and fragmentation of lamellar blocks, whereas cold-drawing of annealed speci- 
mens up to a longitudinal strain of 80% does not induce spatial heterogeneity of their 
deformation; (ii) the elastic modulus increases with the strain rate for non-annealed 
specimens and decreases for annealed samples. 
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1 Introduction 



This paper is concerned with the influence of strain rate on the viscoplastic response of 
isotactic polypropylene (iPP) at isothermal loading with finite strains. The viscoplastic flow 
in semi crystalline polymers, polyethylene [1|, |2|, |3|, f|, |5], ||, |7|, ||, [|, polypropylene JTI], [LI], [12, 
HI, m, m, 0, m, 0, gfj m, H and poly(ethylen e terephthalate) ||, || ||, ||, 0, g| 
|2T| has been a focus of attention in the past decade, which may be explained by the 
importance of the yielding phenomenon for the analysis of fiber- and film-drawing processes. 

Injection- molded isotactic polypropylene is chosen for the experimental analysis for two 
reasons: 

1. Polypropylene is widely used in industrial applications (ranged from oriented films for 
packaging to nonwoven fabrics and reinforcing fibres). 



2. Injection-moded iPP is a semicrystalline polymer with a rather complicated crystalline 
morphology that is strongly affected by annealing JO], [12], |31], [32], [33], |34|, [36| and 

mi 



mechanical loading [10 



21 



Isotactic polypropylene contains three basic crystallographic forms |Pq| : monoclinic a 
crystallites, (pseudo) hexagonal /3 structures, orthorhombic 7 polymorphs, and "smectic" 
mesophase (arrays of chains with a better order in the longitudinal than in transverse chain 
direction). At rapid cooling of the melt (which is typical of the injection-molding process), a 
crystallites and smectic mesophase are mainly developed, whereas (3 and 7 polymorphs are 
observed as minority components |nj. The characteristic size of a spherulites in an injection- 

These structures are formed 



16, 31 



molded specimen is estimated as 100 to 200 /zm 
by lamellae stacks with lamellar thicknesses ranging from 10 to 20 nm [IB, 31] A unique 
feature of the crystalline morphology of iPP is the lamellar cross-hatching: development of 
transverse lamellae oriented in the direction perpendicular to the direction of radial lamellae 
in spherulites [ p5| . 

reveal 



Transmission electron microscopy |32|, |34| and polirized optical microscopy [32 



that annealing of iPP in the sub-melting interval of temperatures results in 

• a decrease in the concentration of transverse lamellae (and disappearance of cross- 
hatching after thermal treatment above 150 °C), 

• an increase in the fraction of ordered a crystallites (where helices in monoclinic unit 
cells are oriented in the same direction). 



WAXS (wide-angle X-ray scattering) measurements [34, 36] demonstrate that annealing 
of isotactic polypropylene in the vicinity of a critical temperature T c w 157 to 159 °C induces 
a second-order phase transition in which the ordered a phase is replaced by a new crystalline 
form with a larger lattice volume and a higher mobility of chains. As molecular mobility 
plays the key role in the viscoplastic behavior of semicrystalline polymers, one can expect 
that this transformation strongly affects the response of iPP. To study this phenomenon in 
detail, two series of uniaxial tensile tests are performed: (i) on specimens not subjected to 
thermal treatment, and (ii) on samples annealed for two days at 160 °C prior to loading. 
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Recent studies on mechanically-induced transformations of the crystalline structure of 
iPP [[[(J [16], [18], Q show that loading a specimen in the sub-yield region of deforma- 



tions results in (i) inter-lamellar separation, (ii) rotation and twist of lamellae, (iii) fine 
slip of lamellar blocks (homogeneous shear of layer-like crystalline structures) |2], f|, , (iv) 
chain slip through the crystals |19|, (v) sliding of tie chains along and their detachment 
from lamellar blocks [|19|| , and (vi) activation of the rigid amorphous fraction (part of the 
amorphous phase located in regions surrounded by radial and transverse lamellae). Strain- 
ing of isotactic polypropylene in the post-yield region of deformations induces (i) coarse 
slip of lamellar blocks (heterogeneous inter-lamellar sliding) and their fragmentation, (ii) 
detachment of chain folds and loops from surfaces of crystal blocks, (iii) disintergation of 
crystalline lamellae into a mosaic structure, and (iv) alignment of broken lamellar blocks 
along the direction of maximal stresses and formation of a fibrillar texture 

To develop a constitutive model for the mechanical response of a semicrystalline poly- 
mer, we apply a method of homogenization. According to this approach |7], a complicated 
micro-structure of isotactic polypropylene is replaced by an equivalent phase whose response 
captures essential features of the mechanical behavior of the semicrystalline polymer. An 
amorphous phase (treated as a network of macromolecules bridged by junctions) is chosen 
as the equivalent phase for the following reasons: 

1. The viscoplastic flow in semicrystalline polymers is "initiated in the amorphous phase 
before transitioning into the crystalline phase" ||. 

2. Sliding of tie chains along and their detachment from lamellae play the key role in the 
yielding phenomenon ||19||. 



A semicrystalline polymer is thought of as an equivalent network of chains bridged by 
junctions (entanglements, physical cross-links on the surfaces of crystallites and lamellar 
blocks): 

• To simplify the analysis, the network is treated as incompressible (observations show 
that the degree of compressibility of polypropylene is rather low: the growth of the 
tensile strain up to 8% causes an increase in the volume strain by less than 1% |37| ). 
The hypothesis regarding incompressibility of semicrystalline polymers at large de- 
formations has been previously applied to the analysis of the mechanical response of 
polyethylene J7], polypropylene [[15], |H| and poly(ethylene terephthalate) [^3], p5| . 

• To exclude from the consideration viscoelastic effects (associated with rearrangement 
of chains), all nodes are treated as permanent (active strands cannot separate from 
junctions and dangling strands cannot merge with the network during the experimental 
time-scale). 

• To describe the viscoplastic flow, the network is assumed to deform non-affinely (junc- 
tions can slide with respect to their reference positions under loading). 

Unlike conventional theories of non-affine networks, two sliding processes are introduced. 
The first reflects motion of nodes in amorphous regions and fine slip of crystallites. This 
process is assumed to occur at any intensity of strains. The other is attributed to coarse slip 
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and fragmentation of lamellar blocks and it takes place in the post-yield region of deforma- 
tions only. The difference between these two types of non-affme deformation is that sliding 
of junctions in amorphous domains does not produce dissipation of energy, whereas coarse 
slip and disintegration of lamellae result in a noticeable entropy production. At uniaxial 
tension, the yield point is associated with the elongation ratio, at which fragmentation of 
lamellar blocks starts. 

The objective of this study is two-fold: 

1. To report experimental data in two series of uniaxial tensile tests on injection-molded 
iPP specimens: the samples in the first series of tests were not subjected to thermal 
treatment, whereas those in the other series of experiments were annealed at 160 °C 
prior to loading. 

2. To develop constitutive equations for the viscoplastic behavior of a semi crystalline 
polymer at finite strains and to determine adjustable parameters in the stress-strain 
relations by fitting observations in tensile tests with the strain rates ranging from 
1.5- 1CT 3 to 6.2- 10 -2 s -1 (practically in the entire region of cross-head speeds employed 
in conventional quasi-static tests). 

Our goal is to assess the influence of thermal pre-treatment and strain rate on the viscoplastic 
response of isotactic polypropylene. 

The exposition is organized as follows. Experimental data in uniaxial tensile tests are 
reported in Section 2. Kinematic relations for sliding of junctions are developed in Section 
3. Kinetic equations for the rate-of-strain tensor (that describes sliding of junctions in amor- 
phous regions and fine slip of lamellar blocks) are introduced in Section 4. Strain energy 
density of a semicrystalline polymer is calculated in Section 5. Constitutive equations for 
an isothermal deformation with finite strains and kinetic relations for coarse slip and frag- 
mentation of lamellar blocks are derived in Section 6 by using the laws of thermodynamics. 
Phenomenological expressions for material functions are proposed in Section 7. The govern- 
ing equations are simplified for uniaxial tension of an incompressible medium in Section 8. 
Adjustable parameters in the stress-strain relations are determined in Section 9 by match- 
ing the experimental data. The effect of strain rate on material parameters is discussed in 
Section 10. Some concluding remarks are formulated in Section 11. 

2 Experimental procedure 

Isotactic polypropylene (Novolen 1100L) was supplied by BASF (Targor). ASTM dumbbell 
specimens were injection molded with length 148 mm, width 10 mm and thickness 3.8 mm. 
Uniaxial tensile tests were performed at room temperature on a testing machine Instron-5568 
equipped with electro-mechanical sensors for the control of longitudinal strains in the active 
zone of samples. The tensile force was measured by a standard load cell. The engineering 
stress, a, was determined as the ratio of the axial force to the cross- sectional area of stress-free 
specimens. The true longitudinal stress, S, was calculated by means of the incompressibility 
condition as 

S = ak, 
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where k is the elongation ratio. 

Two series of mechanical tests were carried out. Every test was performed on a new 
specimen. 

The first series of experiments consisted of 7 tests on specimens that were not subjected to 
thermal treatment with the cross-head speeds 5, 10, 25, 50, 100, 150 and 200 mm/min (which 
correspond to the Hencky strain rates 1.54- 10" 3 , 3.02- 10~ 3 , 7.22- 10~ 3 , 1.49- 10~ 2 , 3.09- 10~ 2 , 
4.56- 10 -2 and 6.23- 10~ 2 s _1 , respectively). The specimens were elongated up to the Hencky 
strain en = 0.3 that noticeably exceeded the strain corresponding to the onset of necking 
at all cross-head speeds. The chosen strain rates ensured nearly isothermal experimental 
conditions, on the one hand, and they allowed the viscoelastic effects to be neglected, on 
the other. The duration of tensile tests before necking of samples varied from 4 s for the 
highest cross-head speed to 180 s for the lowest one. Our recent study demonstrated that 
the amount of stress relaxing during this period did not exceed 10% |38| . 

The other series of experiments was performed on specimens that were annealed in an 
oven for 51 h at 160 °C and cooled by air. To minimize the effect of physical aging, mechanical 
tests were carried out one day after thermal treatment. The series consisted of 7 tests 
that were carried out with the same cross-head speeds as the experiments on non-annealed 
samples. The specimens were stretched up to the maximal Hencky strain, en = 0.6. No 
necking of samples was observed in all experiments, except for the test with the maximal 
strain rate, en = 6.23- 10~ 2 s _1 , where the onset of a weak neck was revealed at the elongation 
ratio k ~ 1.75. 

The true longitudinal stress, E, is plotted versus the elongation ratio, k, in Figures 1 to 
7 for the non-annealed specimens (only the experimental data below the necking points are 
presented) and in Figures 8 to 14 for the samples subjected to thermal pre-treatment. 

A noticeable difference is observed between the engineering stress-engineering strain 
curves for non-annealed and annealed specimens. 

A typical stress-strain diagram for a sample that was not subjected to thermal treatment 
demonstrates an increase in the engineering stress with the engineering strain below some 
strain e y (close to the yield strain, e° = 0.13, provided by the supplier) and a very weak 
decrease (practically a horizontal plateau) in an interval between e y and a strain, e n , at which 
a neck is formed. The presence of a rather wide region of maxima on the engineering stress- 
engineering strain curve does not allow the yield stress and the yield strain for a non-annealed 
specimen to be determined with a high level of accuracy (in agreement with the conventional 
approach in viscoplasticity of solid polymers |39| , the yield point is associated here with the 
point of maximum on the engineering stress-engineering strain diagram). This situation 
is rather typical for semicrystalline polymers that are not subjected to thermal treatment 
prior to mechanical testing: a similar observation was reported in [23] for poly (ethylene 
terephthalate). 

A stress-strain curve for an annealed sample reveals an increase in the engineering stress 
below the yield points, a sharp maximum on the engineering stress-engineering strain di- 
agram, and a monotonical decrease in the engineering stress in the post-yield region of 
deformations. The true stresses, S y , corresponding to the points of maxima on the engi- 
neering stress-engineering strain curves, as well as the appropriate elongation ratios, k y , 
are plotted in Figure 15 versus the Hencky strain rate, en- This figure shows that the true 
yield stress, E y , monotonically increases with en, while the elongation ratio for yielding, k y , 
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is practically independent of the strain rate. It is worth noting that the yield strains for 
non-annealed specimens, e y , exceed the engineering yield strains, e y = k y — 1, for annealed 
specimens by a factor of 2 to 3. 

The experimental data depicted in Figure 15 are fitted by the conventional equation (that 
can be developed within the concept of thermal activation of screw dislocations ||) 

E y = E + E 1 \oge n , (1) 

where log = log 10 and the adjustable parameters S m (m = 0, 1) are determined by the least- 
squares technique. Figure 15 demonstrates fair agreement between the observations and 
their approximation by Eq. (1). It should, however, be noted that the area of applicability 
of the phenomenological relation (1) is rather narrow: it cannot be employed at very low 
strain rates (when the viscoelastic effects become important) and at high rates of straining 
(when rate- induced fracture mechanisms are dominant). 

Our purpose now is to derive constitutive equations for the viscoplastic behavior of a 
semicrystalline polymer that can correctly match the experimental data plotted in Figures 
1 to 14. 



3 Kinematics of sliding 

Denote by r the radius vector of an arbitrary point in the initial state (before external loads 
are applied) and by r(t) its radius vector in the deformed state at time t > 0. Transition from 
the initial state of a network to its actual state is determined by the deformation gradient 

F( ( )^. (2) 

Sliding of junctions in the network with respect to their initial positions is treated as a 
transformation of the reference state, when a point with the initial radius vector ro moves 
to the point with a radius vector r s (t). This transformation is described by the deformation 
gradient 

F - w = -ST- 

It is worth noting that a smooth mapping, r s (t, r ), is determined only locally, which means 
that the new (intermidiate) configuration is not necessary Euclidean . 

A semicrystalline polymer is modelled as an incompressible medium, which implies 
that its deformation from the reference state to the deformed state is isochoric (volume- 
preserving). We suppose that transformation of the initial reference state into the new 
reference state (characterized by the deformation gradient F s ) is volume-preserving as well. 

Transformation of the new reference state into the deformed state is determined by the 
deformation gradient 

F ^> = r§- 

dr s (t) 

The subscript indices "s" and "e" indicate that appropriate deformation gradients describe 
sliding of junctions and elastic deformation (the latter means that the strain energy density 
is a function of F e ). 
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According to the chain rule for differentiation, the tensors F(t), F s (t) and F e (t) are 
connected by the relationship 

F(t)=F e (t)-F B (t), (3) 

where the dot stands for inner product. Formula (3) coincides with the multiplicative presen- 
tation of the deformation gradient proposed in |[23|| , where F e is called the "network stretch" 
tensor and F s is referred to as the "slippage stretch" tensor. 
Differentiation of Eq. (2) with respect to time implies that 

dF _ dv(t) dr(t) 
where v(t) = dv{t)/dt is the velocity vector. Introducing the velocity gradient 

and using Eq. (2), we obtain 

dF 

— (t)=L(t)-F(t). (4) 

Bearing in mind that 

^(t) = -F-\t).^(t).F-\t), 

we find from Eq. (4) that 

— (t) = -F-\t).L(t). (5) 
By analogy with Eqs. (4) and (5), we write 

HF dV^ 1 

-^(t) = L s (t) ■ F s (t), ^-(t) = -F; 1 ^) ■ L s (t), (6) 

where 

L -w = Pi 

dr s (t) 

is the velocity gradient for sliding of junctions, and v s (t) = dr s (t)/dt. 
It follows from Eq. (3) that 

^(t) = ±[F{t) ■ F 7 \t)} = §(t) ■ F;\t) + F(t) ■ ^f(t). 
dt dt l J dt dt 

Substitution of Eqs. (4) and (6) into this equality results in 

dF 

^(t)=L(t)-F c (t)-F c (t)-L s (t). (7) 

The left and right Cauchy-Green tensors for elastic deformation are given by 

B e (t) = F e (t) ■ Fj(t), C.(f) = F c T (t) ■ F c (t), (8) 



7 



where T stands for transpose. We differentiate the second equality in Eq. (8) with respect 
to time, use Eq. (7), and find that 



dt 



(t) = 2FJ (t) • D(t) • F e (t) - L S T (0 • C.(f) - C e (t) • L s (t), 



where 



D(t) = - L(t) +L T (t) 



(9) 
(10) 



is the rate-of-strain tensor. Taking into account that 

^(*) = -c e - 1 (0-^(*)-c e - 1 (0, 

and using Eq. (9), we arrive at the formula 

* (t) = -2F- 1 (t) • D(t) • F-\t)} + c;\t) ■ Lj(t) + L s (t) ■ c;\t). 



dt 



(11) 



The first principal invariant of the right Cauchy-Green tensor, C e (t), reads 

J\{t) — C e (t) : I, 

where I is the unit tensor and the colon stands for convolution. Differentiating this equality 
with respect to time and employing Eqs. (8) and (9), we obtain 



dJi 
~~dt 



where 



(t) = 2 B e (t) : D(t) - C c (t) : D s (t) 



D s (t) = i[L s (t) + Lj(t) 



(12) 



(13) 



is the rate-of-strain tensor for sliding of junctions. For an incompressible medium, the second 
principal invariant of the right Cauchy-Green tensor, C e (t), is given by 

J 2 (t) = C;\t) : I. 

Differentiation of this equality with respect to time and use of Eqs. (8) and (11) imply that 

dJ 2 



dt 



it) = -2 B~\t) : D(f) - C~ e \t) : D„(f) 



(14) 



It follows from Eqs. (12) and (14) that the derivative of an arbitrary smooth function, 
^(^l) -h), of the first two principal invariants of the right Cauchy-Green tensor, C e (t), is 
determined by the formula 



^(MtlMt)) = 2{[$ 1 (t)B c (t)-$ 2 (t)B- 1 (t)] :D(t) 

'Mt)Ce(t)-* 2 (t)C:\t)]:I).(t)\, 



(15) 
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where 

<5>m(t) = -^-{ji(t),J2(t)) (m = l,2). (16) 

To describe strain-hardening of a semicrystalline polymer under active loading in the post- 
yield region of deformations, we introduce the deformation gradient 

f(t)=F s (t)-F s 7 1 (t y ) (17) 

from the reference state at the yielding point to the reference state at time t > t y , where 
t y is the instant when yielding occurs. It follows from Eqs. (6) and (17) that the tensor f 
obeys the linear differential equation 

|(t) = L s (t).f(t), f(t y ) = I. (18) 

By analogy with Eq. (8), we introduce the left and right Cauchy-Green tensors for the 
post-yield transformation of the reference state 

b(t) = f(t)-f T (t), c(t) = f T (t)-f(t). (19) 

Differentiating the second equality in Eq. (19) and using Eqs. (17) and (18), we obtain 

dr 

-(t)=2f T (t)-D s (t).f(t). (20) 



Taking into account that 



f/C " 1 (*) = -c- 1 (*)'^(*)-c- 1 (0, 



dt dt 
we find from Eq. (20) that 

^(t) = -2f- 1 (0-D B (t)-[r 1 (t)] T . (21) 

As any transformation of the reference state is volume preserving, the first two principal 
invariants of the tensor c(t) are given by 

j 1 (t)=c(t):I, m = c- 1 (t):I. (22) 

Differentiating Eqs. (22) with respect to time and using Eqs. (19) to (21), we arrive at the 
formulas 

^(t) = 2b(t) : D.(t), ^(t) = -2h-\t) : B s (t). (23) 

It follows from Eq. (23) that the derivative of an arbitrary smooth function, <f>(ji, jo), of the 
principal invariants of the right Cauchy-Green tensor, c(t), reads 

^{m.m) = 2[&(t)b(f) - Mt)h-\t)} : D s (t), (24) 
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where 

Mt) = |^(ji(t),J 2 (t)) (m = 1,2). (25) 

Equations (19), (22) and (24) are valid for any t > t y . In the sub- yield region of deformations, 
when t < t y , we set formally 

b(t) = c(t) = I, m = j 2 (t) = 3, ^(ji(t),j a (t)) = 0. (26) 

Our aim now is to express the rate-of-strain tensor for sliding of junctions, D s (£), in terms 
of the rate-of-strain tensor for macro-deformation, D(£), and some tensors that characterize 
elastic deformation of a specimen. 



4 Kinetics of fine slip 

A semicrystalline polymer is treated as a strongly heterogeneous network of chains linked 
by junctions. Unlike flI7| , where the spatial inhomogeneity of the network is associated with 
micronecking driven by fragmentation of lamellar blocks, we attribute the heterogeneity of 
the network to an inhomogeneity of interactions between chains in the amorphous phase and 
crystalline lamellae with various lengths and thicknesses. 

The network is thought of as an ensemble of meso-regions with arbitrary shapes and 
sizes. The characteristic length of a MR substantially exceeds the radius of gyration for a 
macromolecule, and it is noticeably less than the characteristic size of a sample. 

Deformation of a specimen induces two sliding processes in the network: 

1. Sliding of junctions between chains with respect to their reference positions in stress- 
free meso-domains. 

2. Sliding of MRs in the ensemble with respect to each other. 
Sliding of nodes in meso-domains of an equivalent network reflects 

• sliding of junctions between chains in the amorphous phase, 

• slippage of tie chains along lamellar surfaces [fL9[j , 

• fine slip of crystalline lamellae (homogeneous shearing of layer-like crystallites and 
small displacements of lamellar blocks with respect to one another) [Q . 

Sliding of MRs with respect to each other describes 

• coarse slip of lamellae (inter-lamellar sliding), 

• fragmentation of lamellae into mosaic blocks linked by tie chains, 

• their alignment along the direction of maximal stresses H. 
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To describe evolution of the microstructure of a semicrystalline polymer driven by simul- 
taneous effects of these two sliding processes, we introduce two rate-of-strain tensors, D f and 
D c , and adopt the conventional assumption that the total rate-of-strain tensor, D s , equals 
the sum of D f and D c , 

(27) 



D s = D f + D c , 
'c" refer to fine and coarse slip, respectively. 



where the subscript indices "f ' and 

We suppose that deformation of a specimen induces sliding of junctions in amorphous 
regions and fine slip of lamellar blocks both in the sub-yield and post-yield regions of de- 
formation. Te rate-of-strain tensor, Df, is assumed to be proportional to the rate-of-strain 
tensor, D, 

D f (t) = a(f)D(t). (28) 

For an incompressible isotropic medium, the coefficient of proportionality, a, depends on 
the principal invariants, J\ and J 2 , of the right Cauchy-Green tensor C c . This coefficient 
vanishes at the zero elastic strain, monotonically increases with elastic deformation, and 
tends to some constant a G [0, 1] (the rate of sliding of junctions for a developed viscoplastic 
flow) at relatively large elastic strains. The inequality a > means that junctions move in 
the direction determined by the macro-strain, whereas the condition a < 1 ensures that the 
rate of the steady flow of nodes does not exceed the rate of macro-strain. 

Assuming the tensor D c (t) to vanish in the sub-yield region of deformations and substi- 
tuting Eqs. (27) and (28) into Eqs. (15) and (26), we find that for any t < t y , 

d *'"-* " -[^(Beit) - a(t)C e (t)) 

-* 2 (t)(B;\t)-a(t)C;\t))]:I>(t), 



dt V 



(Ji(t),j 2 (t)) 



dt 



{Mt)j 2 (tj) = 0. 



(29) 



It follows from Eqs. (15), (24), (27) and (28) that in the post-yield region of deformations, 
when t >t y , the derivatives of the functions $ and read 



d$ ( T / S -r / 



-* 2 (t)(B;\t)- a (t)c;\tj) 



:D(f) 



^{m,m) = 2a(t)[MtMt)- ( f )2 (t)h- 1 (t)} :D(t) 



+2 



MtMt)-Mt)h- L (t) :D c (t), 



(30) 



where the functions $ m and 



1,2) are given by Eqs. (16) and (25). 



5 Strain energy density 

We assume that the macro-strain is transmitted to individual meso-regions by links (crys- 
talline lamellae) that connect MRs in an ensemble. A meso-region is treated as an incom- 
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pressible isotropic medium with the strain energy density (per unit volume) 



w = jJ,W, 



where \i is an average rigidity of a meso-domain and w — w( Ji, J 2 ) is a dimensionless function 
that satisfies the condition 



w(Ji, J 2 ) 



0. 



J 1= 3, J 2 =3 

This equality means that the mechanical energy of a MR vanishes in the reference state. 

It is postulated that the rigidity, fx, is constant in the sub-yield region, and it monotoni- 
cally increases with viscoplastic strains in the post-yield region of deformations, where \x be- 
comes a function of the principal invariants of the right Cauchy-Green tensor c: 11 = fj,(ji, j 2 ). 
This dependence reflects strain-hardening of a semicrystalline polymer induced by texture 
formation. 

We adopt the conventional assumptions that (i) the excluded-volume effect and other 
multi-chain effects are screened for an individual chain in a network by surrounding macro- 
molecules, and (ii) the energy of interaction between chains in a meso-region and between 
meso-domains can be taken into account with the help of the incompressibility condition. 
These hypotheses imply that the strain energy density (per unit volume) of a network, W, 
equals the sum of the mechanical energies of MRs, 

W(t) = M(j\(t), j 2 (0)w(Ji(0, Mtj), (31) 

where M = fiN is the rigidity of the network, and N is the average number of MRs per unit 
volume. 

Differentiating Eq. (31) with respect to time, using Eq. (29), and taking into account that 
the first principal invariant of the rate-of-strain tensor, D, valishes for isochoric deformations, 
we find that in the sub- yield region of deformations, when t < t y , 



dW 
~dT 



(t) = A'(t) : D(f), 



where the prime stands for the deviatoric component of a tensor, 

A(t) = 2M (0)[ Wl (t)(B e (t)-a(t)C e (t)) - w 2 (t)(B; 1 (t) - a(t)C-\t) 



(32) 



(33) 



Bearing in mind that the first principal invariants of the rate-of-strain tensors, D and D c , 
equal zero for a volume-preserving transformation, we find from Eqs. (30) and (31) that in 
the post-yield region of deformations, when t >t y , 



dW 
~dT 



(f)=A'(f):D(f)-A;(f):D c (f), 



(34) 



where 



A(i) = 2{M (t)[w 1 (t)(B e (t) - a(t)C e (t)) - w 2 (t)(B; 1 {t) - aftC; 1 ®) 
+a(t)w (t)\M 1 (t)b(t) - M 2 (t)h-\t) 



A c (t) = 2|M o (0[^i(0CeW-^2(0C e - 1 W]-«;oW[M 1 (t)b(i)-M 2 (i)b- 1 (t)]}. (35) 
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The functions w m (t) and M m (t) (m = 0, 1, 2) in Eqs. (33) and (35) read 



w (t) = w(Ji(t), Mtj), w m (t) = -^-(j 1 (t),J 2 (t)), 

M (t) = M( 3l (t), j2 (t)), Mm (t) = ^(j 1 (t),j 2 (t)). (36) 

Our purpose now is to derive stress-strain relations for a semicrystalline polymer and kinetic 
equations for the evolution of the rate-of-strain tensor, D c (£), in the post-yield region of 
deformations by using the laws of thermodynamics. 

6 Constitutive equations 

For isothermal deformation of an incompressible medium at a reference temperature To, the 
Clausius-Duhem inequality reads |4Cj 



To^(t) = ~(t) + E'(t) : D(t) > 0, (37) 

where XI is the Cauchy stress tensor, and Q is the entropy production per unit volume. 

Substition of Eqs. (32) and (34) into Eq. (37) implies that in the sub-yield region of 
deformations, when t < U 



T ^(f)= [23(f) -A(t)l':D(f), (3* 



dt 

whereas in the post-yield region, when t > t y , 

To^(t) = [S(f) - A(f)]' : D(f) + A^(f) : D c (f). (39) 

The main hypothesis of the present study is that 

• sliding of nodes in the amorphous phase and fine slip of lamellar blocks do not induce 
dissipation of mechanical energy, 

• an increase in the specific entropy of an equivalent network is driven by coarse slip and 
fragmentation of lamellae only. 

This assumption implies that in the sub-yield region of deformations, 

f(*) = 0. (40, 

The entropy production in the post-yield region of deformations is attributed to the coarse 
slip of lamellar blocks characterized by the rate-of-strain tensor D c . Adopting the con- 
ventional formula for the rate of entropy production driven by a viscoplastic flow, we find 
that 

fw = ^ki) D ^ :D <M' < 41 > 
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where 7 is the rate of viscoplastic flow and 



A = (|D : D) 



(42) 



is the intensity of macro-strain rate. 

Substitution of expression (40) into Eq. (38) implies that for an arbitrary loading pro- 
gram, the Clausius-Duhem inequality is satisfied in the sub-yield region of deformations, 
provided that the Cauchy stress tensor, S, reads 

E(f) = -p(t)I + 2M (0)[ Wl (t)(B e (t)-a(t)C e (t)) 

-^(^(B^^-a^C- 1 ^))], (43) 

where pit) is pressure. 

It follows from Eqs. (39) and (41) that for an arbitrary loading program, the Clausius- 
Duhem inequality is satisfied in the post-yield region of deformations provided that the 
Cauchy stress tensor, S, is given by 

E(f) = 



-p(t)I + 2|M (0 [ Wl (t)(B e (t) - a(t)C e (t)) 
+a(t)w (t) [Mi(t)b(t) - M 2 (t)b _1 (t) 



(44) 

and the rate-of-strain tensor for coarse slip of lamellar blocks is determined by the formula 



D c (t) 



2 7 |m (0 



w 1 (t)c e (t)-w 2 (t)c; 1 (t) 

-w (t) M 1 (t)b(t) -M 2 (t)b- 1 (t)]\ , D i (t). 



(45) 



Equation (45) implies that the rate-of-strain tensor for viscoplastic flow in the post-yield 
region, D c , is proportional to the intensity of the rate-of-strain tensor for macro-deformation, 
D. 

The set of constitutive equations for a semicrystalline polymer consists of the stress- 
strain relations (43) and (44) together with the kinetic equations (27), (28) and (45) for fine 
and coarse slip of lamellar blocks. The governing equations are determined by 3 material 
functions: a(t), M(j 1 ,j 2 ) and w( Ji, J 2 ). 

For an arbitrary three-dimensional deformation, the set of governing equations should be 
completed by a relationship between the vorticity tensor for the viscoplastic flow, Q s (t) = 
|[Lj — L s (£)], and the vorticity tensor for macro-deformation, ft(t) = |[L T — L(t)]. We do 
not dwell, however, on such a relation, because this work focuses on uniaxial tension of an 
incompressible medium, when the tensors, f2 and fl s , vanish. 



7 Material functions 

The dimensionless strain energy density, w, is given by 

w = J± — 3. 



(46) 
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Equation (46) describes the mechanical energy of a neo-Hookean medium. The advantages 
of this equation are that it (i) contains no adjustable parameters, and (ii) has a transparent 
physical meaning as the strain energy density of a Gaussian network of flexible chains . 
The following phenomenological relation is proposed for the function a(t): 



a 



1 — cxp 



(47) 



Formula (47) is determined by two adjustable parameters, a and e: the coefficient a is the 
rate of sliding of junctions for a developed viscoplastic flow, and the strain e characterizes 
transition to the steady flow. Equation (47) was successfully employed in our previous study 
on the viscoplastic response of iPP at small strains |j42| . Similar relation was suggested in 
[ f!3|j to describe the time-dependent response of polyethylene melts. 

To describe changes in the elastic modulus, M, induced by texture formation in the 
post-yield region of deformations, we adopt the phenomenological equation 



M(jx,j 2 ) 



-E 



I + tjOi-3) 



(4? 



An advantage of Eq. (48) is that it is determined by two adjustable parameters, E and rj, 
and the effect of viscoplastic flow on the elastic modulus is accounted in a fashion similar to 
that employed in Eqs. (46) and (47). 

As this work focuses on experimental data in uniaxial tensile tests, an explicit expression 
for the yield surface cannot be verified. We accept the von Mises criterion and assume that 
in the sub-yield region of deformations, 



S; < S 



y • 



whereas in the post-yield region, 
where 



Si > S 



y ■ 



S i= (-£': S ,VJ 



(49) 



is the true stress intensity. 

Substituting expressions (46) and (48) into Eqs. (43) to (45) and using Eq. (36), we find 
that in the sub-yield region of deformations, 



£(t) = -p(t)I + E[B c (t)-a(t)C c (t) 
D s (t) = a(t)D(t), 
and in the post-yield region of deformations, 

-p(t)I + e\\i+ rtfj^t) - 3)1 [B e (t) - a{t)C e {t) 



(50) 



E(t) = 

+r ] a(t)[J 1 (t) -3)b(t) 
D„(t) = a(t)D(t)+ 7 EA(t){[l + ^(ji(t)-3)]C c (t) 
-r ? (j 1 (t)-3)b(t)V. 



(51) 
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The stress-strain relations (50) and (51) are determined by 5 adjustable parameters: a, E, 
7, 7] and e. An advantage of this model is that the number of constants to be found by 
fitting experimental data is substantially smaller than that in other constitutive relations in 
finite viscoplasticity of solid polymers E3, SOL §31 §R M, H7L H3 . 



Our purpose now is to simplify the stress-strain relations for uniaxial tension of an 
incompressible medium. 



8 Uniaxial tension 



Points of a bar refer to Cartesian coordinates {X{\ (i = 1,2,3) in the stress- free state, to 
Cartesian coordinates {x{\ in the deformed state, and to Cartesian coordinates in the 
reference state at time t. Uniaxial tension of the incompressible medium is described by the 
formulas 

xx = k(t)X h x 2 = k-*(t)X 2 , x 3 = k-*(t)X 3 , (52) 

where k = k(t) is an elongation ratio. It is assumed that transformation of the reference 
state is determined by the equations similar to Eq. (52), 

& = K (t)X 1 , 6 = «r*(t)*2, 6 = «"'W^3, (53) 
where n(t) is a function to be found. It follows from Eqs. (2), (3), (8), (52) and (53) that 



B r 



-) eiei + ^( e 2e 2 + e 3 e 3 ), 



K 



(54) 



where e« are base vectors of the Cartesian frame {Xi}. By analogy with Eq. (54), we find 
that in the post-yield region of deformations, 



b = c 



K \2 



exe! + — e 2 e 2 + e 3 e 3 

K v 



(55) 



where K y is the value of k at the yield point. It follows from Eqs. (22), (54) and (55) that 



J x (k, k) = (- + 2-, 3i (k, K y ) = (— 

V K 7 k V K V 7 



K \2 

y 



+ 2^. 



(56) 



Substitution of Eqs. (54) and (55) into Eqs. (50) and (51) implies that 

S = Eieiei + S 2 (e 2 e 2 + e 3 e 3 ). 

In the sub- yield region of deformations, the non-zero components of the Cauchy stress tensor, 
Si and E 2 , are given by 

k\2 



51 = -p(k,K)+E\l-a(k,K)](-) 

5 2 = — p(k, k) + E 1 — a(k, k) 



(57) 
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In the post-yield region of deformations, these quantities read 



Si = —p(k, re, K y 



E 



1 - a(k, re)] [l + ??(ji(re, re y ) - 3)] (^) 
+i]a(k, re)(ji(/c, re) - 3) (^-) |, 

S 2 = «, « y ) + #j 1 - cx(k, re) 1 + r](ji(K, Ky) - 3) 

+r/a(k, re)( J\{k, re) — 3) — 



2L 
re 



(58) 



where the function re) is given by Eqs. (47) and (56). 

Excluding the unknown pressure, p, from Eqs. (57) and (58) and the boundary condition, 
E 2 = 0, on the lateral surface of the bar, we find the only component of the Cauchy stress 
tensor, the true longitudinal stress E = Ei. Bearing in mind that for uniaxial tension of an 
incompressible medium, E ; = E, we obtain 



E 
E 



E 
E 



1 — a(k, re) 



/k\ 2 re 
W ~k 



(S < S y ), 



~ir / \ ~i r / A< \ 2 ■ 

1 - a(fc, re) J [l + ^(ji(re, re y ) - 3jJ [(-) - - 



+r)a(k, re)(ji(A;, re) - 3) [(^) 2 - ^] } (E > E y ) 



re \2 re 

Equations (10), (13), (52) and (53) result in the formulas 



(59) 



D 



ldk 
k~dl 



eiei - -(e 2 e 2 + e 3 e 3 ; 



which, together with Eq. (42), imply that 



A 



ldk 

k dt 



1 dn 
re dt 



e^i - -(e 2 e 2 + e 3 e 3 ; 



(60) 



(61) 



It follows from Eqs. (50), (51), (60) and (61) that 



(ire 

~dk 

dn 

Ik 



a{k,n)- (E < E y ), 



-{a(k, re) + T 



-,(*<*, «)- 3 )(£)'-£) 



fc\2 re 



re 



fc 



(S > E y ) 



(62) 



with T = §7f?. 

Given a loading program, /c = the longitudinal stress, S(t), is determined by Eq. 
(59). The elongation ratio, re(i), that characterizes fine and coarse slips of lamellar blocks, 
is found from the nonlinear differential equations (62) with the initial condition re(l) = 1. 
The constitutive equations involve 5 adjustable parameters: 
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1. the elastic modulus E, 

2. the rate of a developed viscoplastic flow of junctions a, 

3. the strain e that characterizes transition to a steady flow of junctions, 

4. the rate of coarse slip of lamellar blocks T, 

5. the coefficient rj that characterizes strain-hardenig in the post-yield region of deforma- 
tions. 

The yield stress, E y , is determined directly from a stress-strain diagram as the true stress 
corresponding to the point of maximum on the engineering stress-engineering strain curve. 
To find the viscoplastic elongation ratio at yielding, « y , the constitutive equations are inte- 
grated from E = to E = E y . 

An important advantage of Eqs. (59) and (62) is that 3 material constants, E, a and e, 
are found by fitting a stress-strain curve below the yield point. Afterwards, the other two 
parameters, T and rj, are determined by matching the stress-strain curve in the post-yield 
region of deformations. 



9 Fitting of observations 

We begin with matching the stress-strain diagrams below the apparent yield point. To find 
the constants E, a and e, we fix some intervals [0, a max ] and [0, £ max ], where the "best-fit" 
parameters a and e are assumed to be located, and divide these intervals into J subintervals 
by the points a® = iAa and = jAe (i, j = 1, . . . , J) with Aa = a max / J and Ae = £ max / J- 
For any pair, {a^ l \e^}, Eqs. (59) and (62) are integrated numerically by the Runge-Kutta 
method with the step Ak = 1.0 • 10~ 5 in the interval between E = and E = E y . Given 
a pair, {a^ l \e^} : the elastic modulus, E, is found by the least-squares method from the 
condition of minimum of the function 

R = 5Z[^cxp(^n) — ^num(^n) > 

where the sum is calculated over all experimental points, k n , in the sub-yield interval, E exp 
is the stress measured in a tensile test, and S num is given by Eq. (59). The "best-fit" 
parameters a and e are determined from the condition of minimum of the function R on the 
set{a« e 0-) (*,j = l,...,J)}. 

The material constants E, a and e that minimize the discrepancies between the exper- 
imental data and the results of numerical analysis are found for any stress-strain curve 
independently. 

First, we approximate the observations for non-annealed specimens. Figures 1 to 7 
demonstrate good quality of fitting the stress-strain curves up to the points of necking 
formation. This result is rather surprising. In terms of the model, it means that for samples 
not subjected to thermal treatment, necking (transition from a homogeneous to a spatially 
heterogeneous deformation of samples) precedes yielding (coarse slip and fragmentation of 
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lamellar blocks). This implies that the presence of a wide plateau on a stress-strain dia- 
gram for a non-annealed specimen near the point of maximum may be attributed not to the 
material yielding (as conventional scenarios for yielding suggest), but to a developed flow of 
junctions in the amorphous phase and fine slip of lamellar blocks in the crystalline phase 
without lamellar fragmentation and texture formation. 

Afterwards, we match the stress-strain curves in the sub- yield region of deformations for 
annealed specimens. Figures 8 to 14 show fair agreement between the experimental data 
and the results of numerical analysis. 

The adjustable parameters E, a and e are plotted versus the rate of Hencky strain, en, 
in Figures 16 to 18. The experimental data are fitted by the functions 

E = E + E x log e H , a = a + a x log e H , log e = e + e 1 log e H , (63) 

where the coefficients E m , a m and e m (m = 0, 1) are found by the least-square technique. 
Figures 16 to 18 demonstrate acceptable agreement between the observations and their 
approximations by phenomenological relations (63). It is worth noting rather large scatter 
of the experimental data for annnealed samples compared to that for specimens not subjected 
to thermal pre-treatment. These discrepancies may be explained by a substantial decrease 
in the interval of elongation ratios where the fitting procedure is performed (the maximal 
engineering strain for the interval where the stress-strain curves are approximated is reduced 
by a factor of three: from 0.15 to 0.05). 

To find the quantities T and 77, we approximate the stress-strain curves for annealed 
specimens in the post-yield region of deformation. For any set of experimental data, we use 
the parameters E, a and e found by fitting an appropriate stress-strain curve in the sub- 
yield region of deformations. To approximate a stress-strain curve above the yield point, we 
apply an algorithm similar to that used to match the observations in the sub-yield region of 
deformations. We fix some intervals [0, r max ] and [0, ?7 max ], where the "best-fit" parameters 
T and rj are assumed to be located, and divide these intervals into J subintervals by the 
points rW = iAT and rj^ = jAq = 1, . . . , J) with Ar = r max /J and Ar/ = i] milx /J. 
Given a pair, {T^\r]^}, Eqs. (59) and (62) are integrated numerically by the Runge-Kutta 
method with the step AA; = 1.0 • 10~ 5 . The "best-fit" parameters T and rj are found from 
the condition of minimum of the function R on the set {T^\r]^ — J)}. 

The dimensionless parameters, T and 77, are plotted versus the Hencky strain rate, en, in 
Figures 19 and 20. The observations for the rate of coarse slip of lamellar blocks show that 
T is constant, T = T . The experimental data for the parameter 77 are approximated by the 
function 

V = Vo + r?iloge H , (64) 

where the coefficients i] m (m = 0, 1) are found by the least-squares algorithm. Figure 20 
demonstrates that Eq. (64) ensures acceptable quality of matching the observations. 

10 Discussion 

Figures 1 to 14 demonstrate fair agreement between the experimental data in uniaxial tensile 
tests and the results of numerical simulation, which means that the constitutive equations 
can be successfully applied to fit observations for isotactic polypropylene (non-annealed 
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as well as subjected to thermal pre-treatment) at various strain rates and elongation ratios. 
This conclusion contradicts the conventional standpoint |41] that (i) the concept of Gaussian 



networks of flexible chains (which implies Eq. (46) for the mechanical energy per unit volume) 
fails to correctly approximate stress-strain curves for solid polymers at finite strains and (ii) 
more sophisticated expressions are necessary for the strain energy density w (based on either 
the slip-link theory flT5| , |T7| , |2"5| or a multi-chain version of the theory of polymeric networks 



with finite extensibility of chains [|7|, |30], [44]]). Within the model proposed, an acceptable 



quality of matching observations is reached based on the hypotheses about (i) non-affine 
motion of nodes in the network and (ii) the growth of the elastic modulus driven by texture 
formation. 

Figure 16 reveals that the elastic modulus, E, monotonically increases with the rate of 
stretching for non-annealed specimens and decreases for annealed samples. The following 
explanations may be provided for this finding. 

As the present study concentrates on the viscoplastic behavior of semicrystalline poly- 
mers, we disregard the viscoelastic phenomena associated with rearrangement of chains (sep- 
aration of active strands from junctions and merging of dangling strands with the network) 
and treat a polymer as an equivalent network of chains bridged by permanent junctions. 
Given a strain rate, the simplest way to account for these effects is to distinguish between 
two groups of meso-regions in an ensemble: 

• the characteristic time for rearrangement of chains in MRs belonging to the first group 
is substentially smaller than the characteristic time of loading (which means that the 
contribution of these meso-domains into the strain energy of a polymer is negligible) , 

• the characteristic time for rearrangement of chains in MRs of the other group exceeds 
the characteristic time for loading (which implies that rearrangement of chains in 
these meso-domains may be disregarded and they may be considered as permanent 
networks) . 

According to this division of MRs into two groups, the number of meso-regions with perma- 
nent junctions per unit volume, N, in the formula for the elastic modulus, E, grows with the 
Hencky strain rate en (because an increase in en results in a decrease in the characteristic 
time of loading, and, as a consequence, a decrease in the content of MRs where stresses in 
chains relax due to the rearrangement process). This conclusion is in agrement with the 
experimental data presented in Figure 16 for non-annealed specimens (curve 1). When the 
cross-head speed decreases from 200 to 5 mm/min (which corresponds to an increase in the 
loading time from 4 to 180 s), the elastic modulus, E, is reduced by about 15%, which 
is quite comparable with the intensity of decay in the longitudinal stress during 180 s in 



uniaxial tensile relaxation tests 38 



The substantial decrease in the elastic modulus of annealed specimens with the Hencky 
strain rate (curve 2 in Figure 16) may be associated with alteration of the morphological 
structure of iPP at thermal treatment. 

Stretching of a semicrystalline polymer in the sub-yield region of deformations activates 
sliding of junctions in amorphous domains between lamellae. Displacements of junctions in- 
duce separation of tie chains from crystallites and their slippage along lamellar blocks. When 
the tie chains are rather long, they have enough time to attach to lamellae in new cites, which 
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implies that the total number of links between amorphous regions and crystalline lamellae 
(than transmit the macro-strain to amorphous domains) remains practically constant. The 
latter means that amorphous meso-domains do not separate from the ensemble at stretching 
of a non-annealed specimen. 

The situation changes dramatically when a semicrystalline polymer is subjected to ther- 
mal pre-treatment. Annealing of iPP for 51 h at 160 °C results in a noticeable increase in 
the degree of crystallinity. Crystallization of polymeric chains in amorphous regions occurs 
strongly inhomogeneously: long tie strands (that have high molecular mobility) are mainly 
crystallized, whereas majority of short tie chains remain in the amorphous state. This implies 
that the average length of a tie strand substantially decreases at thermal treatment. Sliding 
of junctions in amorphous domains caused by straining of a specimen induces separation of 
tie chains from crystallites and their slippage along lamellar surfaces. Unlike non-annealed 
specimens, where tie chains are rather long (which implies that their mobility is high and the 
time necessary for their reconnection is small), for annealed samples, the rate of attachment 
of short tie chains (with low molecular mobility) to new cites on the surfaces of crystallites 
becomes comparable with the strain rate. This means that an increase in the Hencky strain 
rate, en, results in a strong decrease in the concentration of links that transmit the macro- 
strain to amorphous domains. As a consequence, individual MRs separate from the ensemble 
(due to the breakage of links between amorphous meso-domains and surrounding regions), 
which implies that the number of "active" meso-regions per unit volume, N, decreases with 
en- As the elastic modulus, E, is proportional to N, the growth of the strain rate induces 
a decrease the elastic modulus of annealed specimens, in accord with the experimental data 
depicted in Figure 16 (curve 2). 

It is worth noting that at small strain rates (when detachment of tie chains from lamellar 
blocks is not pronounced) the elastic modulus, E, of samples subjected to thermal treatment 
substantially (about by twice) exceeds that of non-annealed specimens (due to the growth 
in the degree of crystallinity). With an increase in the Hencky strain rate, the difference 
between the moduli of annealed and non-annealed samples is reduced, and at relatively high 
strain rates (above e c ~ 5 • 1CT 2 s _1 ), the elastic modulus of non-annealed specimens exceeds 
that of annealed ones. 

Figure 17 shows that the dimensionless rate of a developed flow of junctions, a, mono- 
tonically increases with the strain rate. This phenomenon may be ascribed to breakage of tie 
chains that bridge amorphous regions with lamellar blocks. A decrease in the concentration 
of these chains implies that the mobility of junctions in amorphous MRs grows, which is 
observed as an increase in a with en- At small strain rates (below the threshold value, e c ), 
the rate of a steady viscoplastic flow, a, for non-annealed specimens exceeds that for samples 
subjected to thermal treatment (this observation is ascribed to the secondary crystallization 
at annealing that reduces molecular mobility in meso-regions), whereas at high strain rates 
(that exceed e c ), the inverse picture is revealed (the rate of a developed flow of junctions for 
annealed specimens reaches its ultimate value, a = 1). In agreement with our explanation 
for the decrease in the elastic modulus, E, of annealed specimens with en, this result is 
attributed to breakage of tie chains (which have no time enough for their reconnection). A 
decrease in the content of tie chains implies a strong increase in molecular mobility in amor- 
phous meso-regions, which is reflected by the model as the growth of the rate of viscoplastic 
flow a. 
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According to Figure 18, the strain, e, that characterizes transition to a developed flow of 
junctions, increases with the rate of straining both for non-annealed and annealed specimens. 
It is worth noting that the increase in e is more pronounced for annealed samples than for 
non-annealed ones (however, in the former case, the scatter of experimental data is rather 
large) . 

Figure 19 demonstrates that the rate, T, of coarse slip and fragmentation of lamellar 
blocks in the post-yield region of deformations is independent of the strain rate. This results 
appears to be quite natural: it means that the energy necessary for disintegration of lamellar 
blocks exceeds substentially the elastic energy stored in the semicrystalline polymer, which 
implies that the rate of coarse slip is not affected noticeably by mechanical factors. These 
factors do affect, however, the process of lamellar fragmentation, because, in accord with 
Eqs. (59) and (62), the derivative, dn/dk, is proportional to the longitudinal stress, E. 

Figure 20 reveals that the dimensionless parameter rj (that characterizes strain-hardening 
of a semicrystalline polymer) monotonocally increases with the Hencky strain rate, en- This 
assertion appears to be quite natural. Figure 16 shows that the content of tie chains in 
annealed specimens decreases with the strain rate. As a consequence of this reduction in the 
concentration of tie chains, broken lamellar blocks become more mobile, and they are oriented 
easier along the direction of the longitudinal stress. The strain-rate induced intensification of 
the alignment process for disintegrated lamellar blocks implies a more pronounced increase in 
the stress, E, with the elongation ratio, k, in the post-yield region of deformations observed 
in Figure 20. 

According to Eq. (48) and Figure 20, the growth of the elastic modulus, E, driven by 
strain-hardening of iPP (in tensile tests with the maximal Hencky strain €h = 0.6) is of 
order of few per cent. It is worth noting that the same (qualitatively) level of orientation of 
lamellar blocks was recently observed in SAXS (small angle X-ray scattering) measurements 
on another semicrystalline polymer (polyamide PAll) at uniaxial tension up to the necking 
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11 Concluding remarks 

Two series of uniaxial tensile tests have been performed on isotactic polypropylene at room 
temperature. The first series of experiments is carried out on injection-molded specimens 
non-subjected to thermal treatment. In the other series of tests, the samples were annealed 
for 51 h at the temperature T = 160 °C prior to loading. Each series consists of 7 experiments 
with the cross-hear speeds ranging from 5 to 200 mm/min (which cover practically the entire 
range of cross- head speeds employed in conventional quasi-static tensile tests). Stretching of 
a non-annealed specimen is performed up to the onset of necking. Deformation of a sample 
subjected to thermal treatment is carried out up to the maximal Hencky strain en = 0.6. 

Constitutive equations have been derived for the isothermal viscoplastic behavior of a 
semicrystalline polymer at finite strains. A polymer is treated as an equivalent heterogeneous 
network of chains bridged by permanent junctions. The network is thought of as an ensemble 
of meso-regions linked with each other. Under loading, junctions between chains in MRs 
move with respect to their reference positions (this process reflects sliding of nodes in the 
amorphous phase and fine slip of lamellar blocks). At relatively large elongations, the sliding 
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process is accompanied by displacement of meso-regions with respect to each other (which 
reflects coarse slip and fragmentation of lamellar blocks). 

Unlike the conventional definition of the yield point of a solid polymer at uniaxial tension 
as the point of maximum on the engineering stress-engineering strain diagram, we define the 
yield strain as the strain at which coarse slip and fragmentation of lamellar blocks starts. 
This implies that sliding of junctions in MRs with respect to their initial positions occurs at 
any elongation ratio, whereas sliding of meso-domains takes place in the post-yield region of 
deformations only. 

Stress-strain relations and kinetic equations for coarse slip of lamellar blocks are de- 
veloped by using the laws of thermodynamics. These equations are simplified for uniaxial 
tension of an incompressible medium. The governing equations are determined by 5 ad- 
justable parameters that are found by fitting observations. Fair agreement is demonstrated 
between the experimental data and the results of numerical simulation both for non-annealed 
and annealed specimens. 

The main result of the study is that the stress-strain curves (up to the necking points) 
for specimens not subjected to thermal treatment can be fairly well approximated by the 
constitutive relations where coarse slip and fragmentation of lamellar blocks are not taken 
into account. In terms of the model, this assertion means that necking of non-annealed 
samples occurs in the sub-yield region of deformations and precedes their yielding. On 
the contrary, annealed specimens demonstrate pronounced yield points at relatively small 
elongation ratios, far below the strains corresponding to the onset of necking. 

The following conclusions have been drawn: 

1. The elastic modulus, E, monotonically grows with the Hencky strain rate, en, for non- 
annealed specimens and monotonically decreases for annealed ones. The difference in 
the mechanical response of these samples is attributed to a dramatic decrease in the 
average length of tie strands at thermal treatment. 

2. The rate of a developed viscoplastic flow a (that reflects sliding of junctions in amor- 
phous regions and fine slip of crystalline blocks) and the strain e (that characterizes 
transition to the steady flow of junctions) increase with the strain rate. These obser- 
vations are associated with mechanically-induced breakage of tie chains that results in 
the growth of molecular mobility in amorphous meso-regions. 

3. The rate, T, of coarse slip and fragmentation of lamellar blocks in the post-yield region 
of deformations is independent of mechanical factors. 

4. The dimensionless parameter rj (that characterizes strain-hardening of a semicrystalline 
polymer above the yield point) grows with the Hencky strain rate. This observation 
is explained by enhancement of breakage of tie chains with an increase in the rate 
of straining, and, as a consequence, by an acceleration of alignment of disintegrated 
lamellar blocks along the direction of longitudinal stress. 
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List of figures 



Figure 1: The true stress £ MPa versus the elongation ratio A; in a tensile test with the 
cross-head speed 5 mm/min. Circles: experimental data for a non-annealed specimen. Vector 
indicates the beginning of necking. Solid line: results of numerical simulation 

Figure 2: The true stress £ MPa versus the elongation ratio A; in a tensile test with the 
cross-head speed 10 mm/min. Circles: experimental data for a non-annealed specimen. 
Vector indicates the beginning of necking. Solid line: results of numerical simulation 

Figure 3: The true stress £ MPa versus the elongation ratio A; in a tensile test with the 
cross-head speed 25 mm/min. Circles: experimental data for a non-annealed specimen. 
Vector indicates the beginning of necking. Solid line: results of numerical simulation 

Figure 4: The true stress £ MPa versus the elongation ratio A; in a tensile test with the 
cross-head speed 50 mm/min. Circles: experimental data for a non-annealed specimen. 
Vector indicates the beginning of necking. Solid line: results of numerical simulation 

Figure 5: The true stress £ MPa versus the elongation ratio A; in a tensile test with the 
cross-head speed 100 mm/min. Circles: experimental data for a non-annealed specimen. 
Solid line: results of numerical simulation 

Figure 6: The true stress £ MPa versus the elongation ratio A; in a tensile test with the 
cross-head speed 150 mm/min. Circles: experimental data for a non-annealed specimen. 
Vector indicates the beginning of necking. Solid line: results of numerical simulation 

Figure 7: The true stress £ MPa versus the elongation ratio A; in a tensile test with the 
cross-head speed 200 mm/min. Circles: experimental data for a non-annealed specimen. 
Vector indicates the beginning of necking. Solid line: results of numerical simulation 

Figure 8: The true stress £ MPa versus the elongation ratio A; in a tensile test with the 
cross-head speed 5 mm/min. Circles: experimental data for an annealed specimen. Solid 
line: results of numerical simulation 

Figure 9: The true stress £ MPa versus the elongation ratio A; in a tensile test with the 
cross-head speed 10 mm/min. Circles: experimental data for an annealed specimen. Solid 
line: results of numerical simulation 

Figure 10: The true stress £ MPa versus the elongation ratio A; in a tensile test with the 
cross-head speed 25 mm/min. Circles: experimental data for an annealed specimen. Solid 
line: results of numerical simulation 

Figure 11: The true stress £ MPa versus the elongation ratio A; in a tensile test with the 
cross-head speed 50 mm/min. Circles: experimental data for an annealed specimen. Solid 
line: results of numerical simulation 

Figure 12: The true stress £ MPa versus the elongation ratio A; in a tensile test with the 
cross-head speed 100 mm/min. Circles: experimental data for an annealed specimen. Solid 
line: results of numerical simulation 

Figure 13: The true stress £ MPa versus the elongation ratio A; in a tensile test with the 
cross-head speed 150 mm/min. Circles: experimental data for an annealed specimen. Solid 
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line: results of numerical simulation 

Figure 14: The true stress S MPa versus the elongation ratio k in a tensile test with the 
cross-head speed 200 mm/min. Circles: experimental data for an annealed specimen. Vector 
indicates the beginning of necking. Solid line: results of numerical simulation 

Figure 15: The true yield stress S y MPa (unfilled circles) and the elongation ratio for 
yielding k y (filled circles) versus the strain rate iu s _1 . Symbols: treatment of observations. 
Solid lines: approximation of the experimental data by Eq. (1) with So = 39.27 and Si = 
2.32 (curve 1) and by the constant k y = 1.0473 (curve 2) 

Figure 16: The elastic modulus E GPa versus the strain rate en s _1 . Symbols: treatment 
of observations. Unfilled circles: non-annealed specimens; filled circles: annealed specimens. 
Solid lines: approximation of the experimental data by Eq. (63). Curve 1: E = 0.81, 
E 1 = 0.07; curve 2: E = 0.43, E x = -0.22 

Figure 17: The dimensionless rate a of a developed flow of junctions versus the strain rate 
en s -1 . Symbols: treatment of observations. Unfilled circles: non-annealed specimens; filled 
circles: annealed specimens. Solid lines: approximation of the experimental data by Eq. 
(63). Curve 1: a = 0.93, a x = 0.033; curve 2: a = 0.82, a x = 0.042 

Figure 18: The strain e that characterizes transition to a developed flow of junctions versus 
the strain rate en s -1 . Symbols: treatment of observations. Unfilled circles: non- annealed 
specimens; filled circles: annealed specimens. Solid lines: approximation of the experimental 
data by Eq. (63). Curve 1: e = -1.15, e l = 0.028; curve 2: e = -0.62, e 1 = 0.375 

Figure 19: The dimensionless rate of coarse slip of lamellar blocks T versus the strain rate 
e H s _1 . Circles: treatment of observations for annealed specimens. Solid line: approximation 
of the experimental data by the constant r = 6.07 

Figure 20: The dimensionless parameter rj that characterizes strain-hardening versus the 
strain rate en s" 1 . Circles: treatment of observations for annealed specimens. Solid line: 
approximation of the experimental data by Eq. (64) with r] = 0.0328 and rji = 0.0072 
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